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ABSTRACT 

The initial mass function determines the fraction of stars of different intial mass born 
per stellar generation. In this paper, we test the effects of the integrated galactic initial 
mass function (IGIMF) on the chemical evolution of the solar neighbourhood. The 
IGIMF (Weidner & Kroupa 2005) is computed from the combination of the stellar 
intial mass function (IMF), i.e. the mass function of single star clusters, and the 
embedded cluster mass function, i.e. a power law with index j3. By taking into account 
also the fact that the maximum achievable stellar mass is a function of the total mass 
of the cluster, the IGIMF becomes a time-varying IMF which depends on the star 
formation rate. We applied this formalism to a chemical evolution model for the solar 
neighbourhood and compared the results obtained by assuming three possible values 
for /3 with the results obtained by means of a standard, well-tested, constant IMF. In 
general, a lower absolute value of /3 implies a natter IGIMF, hence a larger number of 
massive stars and larger metal ejection rates. This translates into higher type la and II 
supernova rates, higher mass ejection rates from massive stars and a larger amount of 
gas available for star formation, coupled with lower present-day stellar mass densities. 
Lower values of /? correspond also to higher metallicities and higher [a/Fe] values at 
a given metallicity. We consider a large set of chemical evolution observables and test 
which value of /3 provides the best match to all of these constraints. We also discuss 
the importance of the present day stellar mass function (PDMF) in providing a way 
to disentangle among various assumptions for /3. Our results indicate that the model 
adopting the IGIMF computed with (3 ~ 2 should be considered the best since it 
allows us to reproduce the observed PDMF and to account for most of the chemical 
evolution constraints considered in this work. 

Key words: Galaxy: interstellar medium; Galaxies: evolution; Galaxies: abundances; 
Galaxies: star clusters. 



1 INTRODUCTION 

The initial stellar mass function is one of the major in- 
gredients of chemical evolution models. Moreover, observed 
chemical abundances allow one to put robust constraints 
on both the normalization and the slope of the initial mass 
function (IMF; Chiappini et al. 2000; Romano et al. 2005). 
The environment providing most observational constraints 
for chemical evolution studies is the solar neighbourhood (S. 
N. hereinafter), for which a large set of observables are avail- 
able. These observables include diagrams of abundance ra- 
tios versus metallicity, particularly useful when they involve 
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two elements synthesised by stars on different timescales. 
An example is the [a/Fe] vs [Fe/H] diagram, since a ele- 
ments are produced by massive stars on short (< 0.03 Gyr) 
timescales, while type la supernovae (SNe) produce mostly 
Fe on timescales spanning from 0.03 Gyr up to one Hubble 
time (Matteucci 2001). This diagnostic is a strong function 
of the IMF, but depends also on the assumed star forma- 
tion history (Matteucci 2001; Calura et al. 2009). Another 
fundamental constraint is the metallicity distribution of liv- 
ing stars, which provides us with fundamental information 
on the IMF and on the infall history of the studied system. 
Another diagnostic, depending both on the IMF and the 
past star formation history, is the present-day mass func- 
tion, which represents the mass function of living stars ob- 
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served now in the Solar Vicinity (Elmegreen & Scalo 2006). 
Other important observables useful for chemical evolution 
studies include the type la and type II SN rates, as well as 
the surface density of stars and gas, depending on the IMF 
and on the rate at which the gas has been processed into 
stars and remnants in the past, i.e. on the SFR. 

In a previous paper, Recchi et al. (2009) considered a 
star-formation dependent IMF, called the integrated galac- 
tic initial mass function (IGIMF), which originates from the 
combination of the stellar IMF within each star cluster and 
the embedded cluster mass function. ' Within each star clus- 
ter, the IMF can be well approximated by a two-part power- 
law form, £(m) oc m~ Q (e.g. Pflamm-Altenburg, Weidner &: 
Kroupa 2007). Massey &: Hunter (1998) have shown that 
for stellar masses m > a few Mq, a slope similar to the 
Salpeter (1955) index (i.e. a — 2.35) can approximate well 
the IMF in clusters and OB associations for a wide range of 
metallicities. Other studies have shown that the IMF flat- 
tens out below m ~ 0.5 Mq (Kroupa, Tout & Gilmore 1993; 
Chabricr 2003). On the other hand, the embedded cluster 
mass function is well approximated by a single slope power 
law. This implies that small embedded clusters are more nu- 
merous in galaxies and they lock up most of of the stellar 
mass. However, the most massive stars tend to form prefer- 
entially in massive clusters (Weidner & Kroupa 2006). The 
integrated IMF in galaxies, the IGIMF, is a function of the 
galactic star formation rate (SFR) and, as a consequence 
of the embedded cluster mass function, it is steeper than 
the stellar IMF within each single star cluster (Kroupa & 
Weidner 2003; Weidner & Kroupa 2005). 

Recchi et al. (2009) studied the effects of the IGIMF on 
the evolution of the SN rates in galaxies and on the chemical 
evolution of elliptical galaxies, showing how the IGIMF nat- 
urally accounts for the relation between the [a/Fe] and the 
stellar velocity dispersion observed in local elliptical galax- 
ies. In this paper, we consider the effects of the IGIMF on 
the chemical evolution of the solar neighbourhood. As al- 
ready stressed, the advantage of this approach is the avail- 
ability of a large set of observational constraints, useful to 
test the IGIMF and, most importantly, to constrain its main 
parameter, i.e. the index /3 of the power law expressing the 
embedded clusters mass function. 

We will compare the results computed by means of a 
standard IMF, similar to the one by Scalo (1986), success- 
ful in reproducing most of the chemical evolution properties 
of the Solar Neighbourhood, with the results computed by 
means of the IGIMF. This Paper is organized as follows. 
In Section 2 we describe the IMF used in standard chemi- 
cal evolution models and the formalism behind the IGIMF. 
In Section 3 we present a brief description of the chemical 



t Embedded clusters arc stellar clusters that arc partially or 
fully encased in interstellar gas and dust within molecular clouds, 
therefore often visible only in the infrared. It is supposed that all 
(or the large majority of) the stars form originally in embedded 
clusters (Lada & Lada 1991), but then they can loose their cocoon 
of gas because of the feedback of O stars (see e.g. Boily & Kroupa 
2003a, b). We will name therefore hereafter "embedded clusters" 
also the clusters which have lost their envelope, but still retain 
all of their stars, to be consistent with the terminology used in 
the original papers describing the IGIMF theory (e.g. Kroupa & 
Weidner 2003; Weidner & Kroupa 2005). 



evolution model of the Solar Neighbourhood. In Sect. 4 we 
present our results and in Sect. 5 we draw our conclusions. 



2 THE INITIAL MASS FUNCTION 

2.1 The standard initial mass function 

£,std(m) is the initial mass function assumed in the standard 
chemical evolution model used in this paper and is a two- 
slope power law, defined in number as: 



&>td(m) 



0.19 • ra 
0.24 • m~ 



if m < 2 M Q 
if m > 2 Mq, 



(1) 



In the remainder of the paper, we will refer to the IMF 
of Eq. 1 as to the Standard IMF. This equation repre- 
sents a simplified two-slope approximation of the actual 
Scalo (1986) IMF, similarly to what is done in Matteucci 
& Francois (1989). Our basic IMF is assumed to be con- 
stant in space and time and normalized in mass to unity in 
the mass interval 0.1 — 1OOM , i.e.: 



/ 



m£ s td(m)dm = 1. 



(2) 



Various papers have shown that, by assuming this IMF it is 
possible to reproduce a large number of observational con- 
straints for the solar neighbourhood (Chiappini et al. 2001; 
Romano et al. 2005). By means of our chemical evolution 
model, we will present predictions for various observables 
computed by assuming the standard IMF. These predictions 
will be tested against observed quantities and will be com- 
pared to results computed assuming the integrated galactic 
IMF, which is the subject of the following section. 



2.2 The integrated galactic initial mass function 

The IGIMF theory has been described in detail in previous 
papers (Kroupa & Weidner 2003; Weidner & Kroupa 2005; 
Recchi, Calura & Kroupa 2009). Here we briefly summarize 
its main assumptions and features. The IGIMF theory is 
based on the assumption that all the stars in a galaxy form 
in star clusters. Within each embedded cluster, the stellar 
IMF has the canonical form £(m) = kmT a , with a — 1.3 for 
~ 0.1 M Q < m < 0.5 M and a = 2.35 (i.e. the Salpeter 
slope) for 0.5 M < m < m max . The upper mass m max 
depends on the mass of the embedded cluster M cc \ simply 
because small clusters do not have enough mass to produce 
very massive stars. 

On the other hand, star clusters are also apparently dis- 
tributed according to a single-slope power law, £ oc i oc M~f 
(Zhang & Fall 1999; Lada & Lada 2003). In this work we 
have assumed 3 possible values of j3: 1.00, 2.00 and 2.35 . 
By convolving the stellar IMF with the distribution of em- 
bedded clusters we obtain the IGIMF, namely the IMF in- 
tegrated over the whole population of embedded clusters 
forming in a galaxy as a function of the star formation rate 



CiGiMF(m;^(t)) 



,WW) 

£(ra < m max )£eci(M e ci)<iMeci,(3) 



where M cc i jm i n and M ec i,max(^(i)) are the minimum and 
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maximum possible masses of the embedded clusters in a pop- 
ulation of clusters, respectively, and m m ax = Jn max (M cc i). 
For M cc i jm in we take 5 M (the mass of a Taurus- Auriga 
aggregate, which is arguably the smallest star-forming "clus- 
ter" known) . The upper mass of the embedded cluster popu- 
lation depends instead on the SFR and that makes the whole 
IGIMF dependent on tp. The correlation between M cc i jmax 
and SFR has been determined observationally (Larsen & 
Richtler 2000; Weidner et al. 2004) and results from the 
sampling of clusters from the embedded cluster mass func- 
tion given the amount of gas mass being turned into stars 
per unit time (Weidner et al. 2004). 

In Fig. 1, we show the IGIMF as a function of the 
SFR for the three values of j3 considered in this work, com- 
pared to our standard IMF. The IGIMFs are characterized 
by a nearly uniform decline, which follows approximately a 
power law, and a sharp cutoff when m gets close to m ma x- 
Of course, the steepest distribution of embedded cluster (in 
our case the model with /3 = 2.35) produces also the steepest 
IGIMF because this distribution is biased towards embed- 
ded clusters of low mass, therefore the probability of finding 
high mass stars in this cluster population is lower. More- 
over, the dependence of the IGIMF on the SFR is strong for 
SFR < 1 M yr" 1 whereas it is very mild for SFR > 1 M 
yr -1 (sec Rccchi et al. 2009). This is due to the fact that for 
SFRs larger than ~ 1 M yr~ the maximum possible mass 
of the embedded cluster is very high, therefore it is always 
possible to sample massive stars in the whole galaxy up to a 
mass very close to the empirical limit (which is assumed to 
be 150 M Q ; see Weidner & Kroupa 2005). All the IGIMFs 
are normalized in mass to unity as the standard IMF (see 
eq. 2). 

The standard IMF is steeper in the low mass range, 
i.e. for stellar masses m < 0.5 M . These masses do not 
contribute to the chemical enrichment of the Galaxy, since 
they have lifetimes larger than the Hubble time. However, 
their distribution remains unchanged during the evolution 
of the Galaxy and it can be constrained by analysing the 
present-day mass function (see Sect. 4.8). 



3 THE CHEMICAL EVOLUTION MODEL 

The adopted chemical evolution model is calibrated in order 
to reproduce a large set of observational constraints for the 
Milky Way galaxy (Chiappini et al. 2001). The Galactic disc 
is approximated by several independent rings, 2 kpc wide, 
without exchange of matter between them. The Milky Way 
is assumed to form as a result of two main infall episodes. 
During the first episode, the halo and the thick disc are 
formed. During the second episode, a slower infall of external 
gas forms the thin disc with the gas accumulating faster in 
the inner than in the outer region ("inside-out" scenario, 
Matteucci & Frangois 1989). The process of disc formation 
is much longer than the halo and bulge formation, with time 
scales varying from ~ 2 Gyr in the inner disc to ~ 7 Gyr in 
the solar region and up to 20 Gyr in the outermost disc. In 
this paper, we are interested in the effects of a time-variyng 
IMF in the Solar Neighbourhood. For this purpose, we focus 
on a ring located at 8 kpc from the Galactic centre, 2 kpc 
wide. 



For a single-phase gas, the chemical evolution of a given 
chemical element i is computed through the following equa- 
tion: 

dGi(t) 



dt 



-&*(t)Xi{t) 



+ / o~*(t - T m )Q mi (t - r m )£ (m,t)dm 
Jm l 

pM Bm ,0.5 

+A Ia / / /(M)Qmi(t - T m2 )£ (m,i)<7*(i - T m2 )dfidm 

J M Bm J Mmin 

+(1 - A Ia ) / <j*(t - T m )Q mi (t - r m )£ (m,t)dm 
rMu(t) / dGi(t) 

{t - Tm)Qmi{t - T m )i (m,t)dm +( ^ )inf, 



r-M v (t) 



where d(t) = a g (t)Xi(t)/atot is the gas surface mass den- 
sity in the form of an element i normalized to a total surface 
mass density a to t, and G(t) — a g (t)/atot is the total frac- 
tional mass of gas present in the galaxy at the time t. Xi(t) 
is the abundance by mass (or mass fraction) of the element 
i. The quantity <r» (t) is the surface SFR density. In general, 
this quantity is a function of the galactic radius r: 

(Chiappini et al. 1997), where v is the SF efficiency, <r(r, i) 
is the total mass (gas + stars) surface density at a radius r 
and time t (tc a i = 14 Gyr is the present age of the Milky 
Way), a(r & ,t) is the total surface mass density in the solar 
region and a g (r,t) is the ISM surface mass density. In this 
paper, we focus on the solar neighbouhod and we assume 
r — r = Skpc. For the gas density exponent k a value of 
1.5 has been assumed by Chiappini et al. (1997) in order to 
ensure a good fit to the observational constraints for a large 
set of local spirals (Kennicutt 1998). The efficiency of SF is 
set to v — 1 Gyr -1 , and becomes zero when the gas surface 
density drops below a certain critical threshold. For the SF, 
we adopt a threshold gas density a t h ~ 7M pc -2 in the disc 
as suggested by Kennicutt (1989). 

The second term in Eq. 4 is the rate at which each 
element is restored into the ISM by single stars with masses 
in the range Ml - Ms m , where Ml is the minimum mass 
contributing, at a given time t, to chemical enrichment (the 
minimum is O.8M ) and Ms m is the minimum mass allowed 
for binary systems giving rise to type la SN (3M , Matteucci 
& Greggio 1986). The quantity 



£ (m,t) = m£(m,t) 



(6) 



is the initial mass function in mass and in the standard 
case is given by £ std = m^ at d, described in Sect. 2.1 and 
is constant in time, otherwise it is a function of the initial 
stellar mass m and of the time t and is computed with the 
method described in Sect. 2.2. The quantities Q m i(t — r m ) 
(where r m is the lifetime of a star of mass m) contain all the 
information about stellar nucleosynthesis for elements ei- 
ther produced or destroyed inside stars or both (Talbot and 
Arnett 1971; Matteucci 2001). The third term represents 
the enrichment due to binaries which become type la SNe, 
i.e. all the binary systems with total mass between Ms m 
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log (m/M Q ) 

Figure 1. IGIMFs for different cluster mass functions distributions. Upper panel: j3 = 1; central panel: /3 = 2; lower panel: /3 = 2.35. In 
each panel we have considered 7 possible values of SFRs, ranging from 10 -4 Mq j/r -1 (lowermost solid lines) to 100 Mq yr~ 1 (uppermost 
solid lines), equally spaced in logarithm. The thick dashed line is the standard IMF used in this work (see Sect. 2.1). The IGIMFs and 
the standard IMF are all normalized through the equation J m£(m)dm = Mtot = IMq. 



and Mb m = 16Mq. For the type la SN progenitor model, 
the Single Degenerate (SD) scenario is assumed, where a C- 
O white dwarf explodes by C-defiagration mechanism after 
having reached the Chandrasekhar mass (I.44M0), owing to 
progressive mass accretion from a non-degenerate compan- 
ion (Whelan & Iben 1973). This model is still one of the best 
to reproduce the majority of the properties of local galaxies 
(Matteucci et al. 2006; Calura & Matteucci 2006; Matteucci 
et al. 2009). The parameter Ai a represents the unknown 
fraction of binary systems with the specific characteristics 
to become type la SNe in the range 3-16 Mq and is fixed by 
reproducing the observed present time SN la rate (Calura 
& Matteucci 2006). In this third term, both quantities <r, 
and Qmi refer to the time t — t m2 , where t m2 indicates the 
lifetime of the secondary star of the binary system, which 
regulates the explosion timescale. /j, = M 2 /Mb is the ratio 
between the mass of the secondary component M2 and the 
total mass of the binary system Mb, whereas f(fj,) is the dis- 
tribution function of this ratio. Statistical studies indicate 
that mass ratios close to 0.5 are preferred, so the formula: 

/( M ) = 2 1+ ^l + 7 )^ (7) 

is commonly adopted, with 7 = 2 (Matteucci & Recchi 
2001). fimin is the minimum mass fraction contributing to 
the SNIa rate at the time t, and is given by 

/ M 2 (t) M 2 -0.5M B \ , q . 
Mmin = max I __ _ I . ( 8) 



The fourth term represents the enrichment due to stars in 
the mass range Mg m - Mb m which are either single, or, 
if in binaries, do not produce a SN la event. In this mass 
range, all the stars with masses m > 8Mq will explode as 
type II SNe, which in our picture are assumed to originate 
from core collapse of single massive stars. The fifth term 
of Eq. 4 represents the enrichment of stars more massive 
than Mb m , all of which explode as core collapse (i.e mostly 
type II, see Calura & matteucci 2006) SNe. The upper mass 
limit contributing to chemical enrichment is Muu\ and is a 
function of time through the IMF, £ (m, t). We assume that 
the maximum value for Mu(t) is 100 Mq. The assumption 
of a maximum stellar mass 150 Mq would have a negligible 
effect on the results presented in this paper. Finally, the last 
term accounts for the infall of external gas. 

It is worth stressing that equation 4 differs substan- 
tially from the classic formalism used in most of the previous 
chemical evolution models (e.g. Matteucci & Greggio 1986; 
Chiappini et al. 2001), which use an IMF constant in time, 
whereas in this case the IMF is a function of the SFR, which 
is a function of cosmic time. Later on, we will see how the 
time variations of the IGIMF depend on the star formation 
history for various values of the parameter j3. 

The rate at which the thin disc is formed out of external 
matter is 

(^^) i „ f = B(i?) e - (t - t — )/TD , (9) 

where t max is the time of maximum gas accretion onto the 
disc, corresponding to the end of the halo-thick disc phase 
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and it is equal to 1 Gyr. The quantity td is the timescale 
for mass accretion onto the thin disc component. Following 
Romano et al. (2000) and Chiappini et al. (2001), we assume 
that td increases with increasing Galactic radius 



td(R) = 1.033 x R- 1.267. 



(10) 



The constant B(Rq) is fixed in order to reprodue the 
present-day total surface mass density (stars + gas) in the 
solar neighbourhood. Concerning the nucleosynthesis pre- 
scriptions, we assume the yields of Van den Hoeck & Groe- 
newegen (1997) for low and intermediate mass stars, the 
yields of Iwamoto et al. (1999) for type la SNe and those 
of Frangois et al. (2004) for massive stars. In Table 1, we 
show the standard assumptions for the main parameters of 
our chemical evolution model. 



3.1 Observational data set 

In Table 2, we show the solar neighbourhood observables 
used in this paper, with their values and references. 
Concerning the observed type la and type II SN rates, as 
observational data we use the values of Capellaro (1996), 
expressed in century^ 1 and valid for the whole MW disc. To 
compare them with the predictions, we divide the observed 
values by the area of the MW disc, which, assuming a radius 
of ~ 15 kpc, is roughly Sdisc ~ 10 9 pc 2 . In this way, for 
the type la and type II SN rates we obtain values of 0.003 
pc~ 2 Gyr _1 and 0.012 pc _2 Gyr _1 , respectively. 

Concerning the stellar surface density of visible stars 
and remnants, as observational value we use an indirect esti- 
mate. In a recent paper, Weber & de Boer (2009) discuss the 
uncertainty of the total local visible surface density. They 
provide a value of 48 ± 9 Mq /pc 2 . This value is in agreement 
with another recent estimate of Holmberg & Flynn (2004), 
which found 53 M & /pc 2 in visible matter, and is consis- 
tent with the value we achieve with our standard model, 50 
M e /pc 2 . 

The local gas surface density is between 7 and 14 
Mq/pc 2 (Kulkarni & Heyles 1987, Dame 1993, Oiling & 
Merrifield 2001). A reasonable value based on this data is 
10.5 ± 3.5 Mq/pc 2 . This value is also consistent with our 
estimate from the standard model. By combining these two 
quantities, i.e. by subtracting the gas density from the to- 
tal density and by combining the errors, for the local mass 
density in stars and remnants we obtain 37.5 ± 10M Q /pc 2 . 
This estimate is compatible with previous values based on 
the combination of the observed stellar mass density and 
the one in stellar remnants (Gilmore et al. 1989, Mera et al. 
1998). 



3.2 IGIMF and Star Formation Rate 

In Section 2, we have seen that the IGIMF is a function of 
the total SFR ip{t), expressed in M Q /j/r. In particular, the 
upper mass limit of the IGIMF is sensitive to the value of 
the star formation rate. 

Our chemical evolution code is designed to produce all phys- 
ical quantities as surface mass densities, so the direct output 
of our code is a SFR surface density (the quantity a, defined 
in Sect. 3). To convert the SFR density into the appropri- 
ate units, we perform the following assumptions. The solar 



neighbourhood is described by a 2 kpc-wide ring located at 
the galactocentric distance of 8 kpc from the centre. The 
area of this ring is So ~ 10 s pc 2 . The SFR ip(t) can be 
calculated from the SFR surface density a* as 



ip(t) — <t* • Sq • 10" 



(11) 



In this way, we obtain for our standard model a present 
value of 0.26Mq /yr . If we consider that the local observed 
SFR density is 3.5M Q pc~ 2 yr~ 1 and we convert this value 
into the same units, we obtain O.35M /yr, consistent with 
the above estimate. 



4 RESULTS 

Our aim is to test the effects of the IGIMF on the chemi- 
cal evolution of the S.N., taking into account a large set of 
available observational constraints. First, we will consider 
the effects of the SFR-dependent IGIMF on the predicted 
physical properties of the solar neighbourhood. Then, we 
will discuss all of the observables and the main parameters 
which can be tested by our analysis. Finally, we will see how 
a fine-tuning of the parameters considered in this work may 
allow us to derive constraints on the local IMF. 



4.1 The fitness test 

In order to quantitatively compare our results with the ob- 
servables considered in this paper, we adopt the formula 



fitness = 



1 



1+5 



w(i)[obs(i) — theo(i)] 2 
max {[obs(i)] 2 , [theo(i)] 2 } 



(12) 



where, for the i-th value of each considered parameter, obs 
(i) and theo (i) are the observed values and the predictions 
of the model, respectively (see also Ruzicka et al. 2007; Theis 
& Kohle 2001). We have also introduced a weight w(i) for 
each observable i in order to give each set of observables 
the same statistical weight. To be more precise: we have 
grouped the observables in 5 groups: (1) physical quantities 
(SFR, SNR//, SNR/ a , E», S gas ); (2) solar abundances; (3) 
average [a/Fe] ([O/Fe], [Si/Fe]) ratios for each bin of [Fe/H]; 
(4) dN/d[Fe/H] (SMD) in each bin of [Fe/H]; (5) present-day 
mass function in each bin of mass. We have chosen a weight 
such that w(i) multiplied by the members of each group 
gives always the same number. This means that in the fitness 
calculation, the same weight is given to all the 5 groups of 
observables. Following this approach, each single observable 
in the physical quantities group has the largest weight, and 
this occurs because the physical quanties group is the one 
having the smallest number of members (5). On the other 
hand, each single member of the [a/Fe] group, counting the 
highest member number (28), has the lowest one, but, as 
already stressed, the [a/Fe] group as a whole has the same 
weight of the physical quantities group. Of course, the closer 
fitness is to 1, the better the model is in reproducing the 
observations. In Fig. 4, we show the "fitness" quantity as 
a function of /? for all the models considered in this paper. 
The results obtained for each value of P will be discussed 
separately. 



© — RAS, MNRAS 000, 1-22 



6 F. Calura, S. Recchi, F. Matteucci, P. Kroupa 



Table 1. Parameters for the standard model of the solar neighbourhood. 



Parameter 


Adopted value 




SF efficiency 






v (Gyr -1 ) 


1 




Infall timescalc for the 






thin disc (Gyr) 


7 




Gas density threshold 






for SF 






(M /pc 2 ) 


7 




Fraction Aj a of binary systems 






originating type la SNe 


0.04 




Table 2. Solar neighbourhood observables used in this paper 


and references. 


Observable 


value 


Reference 


SFR Surface density 


3.5 ±1.5 M Q pc- 2 Gyr- 1 


Rana (1991) 


type la SNR 


0.003 ±0.002 pc^Gyr- 1 


Cappellaro (1996) 


type II SNR 


0.012 ±0.008 pc^Gyr- 1 


Cappellaro (1996) 


Gas surface density 


10.5 ±3.5 M Q pc- 2 


Kulkarni & Heiles (1987) 






Dame (1993) 






Oiling & Merrifield (2001) 


Stellar surface density 


37.5 ± 10 M Q pc- 2 


Weber & de Boer (2009) 


(visible stars and remnants) 






Stellar abundance ratios 




various authors 


Stellar Mctallicity distribution 




various authors 


Present-day mass function 




various authors 



4.2 The effects of the IGIMF on the properties of 
the solar neighbourhood 

In Fig. 2 we show the effects of the standard IMF and of the 
IGIMF assuming three different values for /3 on the calcu- 
lated evolution of the star formation rate, of the type la and 
II SN rates and of the stellar and gas mass surface density. 

All the models shown in Fig. 2 are characherised by the 
same value for the star formation efficiency, v = 1 Gyr -1 . 
In this way, it is possible to appreciate the effects of varying 
the IMF, keeping all the other parameters constant. 

The most striking features of the standard model are 
the SF hiatus at 1 Gyr (see Chiappini et al. 1997, 2001) and 
the threshold-dominated SF after 10 Gyr. 

The hiatus in the SF is due to the transition between 
the end of the halo/thick-disc phase and the beginning of 
the thin-disc phase. This effect is confirmed by the relation 
between [Fe/O] and [O/H] observed in local stars (Chiap- 
pini et al. 1997; Gratton et al. 2000) and can be naturally 
reproduced once a SF threshold is adopted (Chiappini et al. 
2001). 

Moreover, in the standard model the adoption of a 
SF threshold causes numerous oscillations in the calculated 
SFH, which is also reflected in the SNRs. 

Within the first 2 Gyr of evolution, the SFRs of the 
standard model and the ones calculated for the models 



adopting the IGIMF are all very similar. The SF hiatus is 
visible also in the models with the IGIMF, and this indi- 
cates that the threshold effect is important in all the four 
cases studied in Fig. 2. For times larger than 2 Gyr, the 
SFRs for models computed with the IGIMF are not very 
sensitive to the adoption of the SF threshold, with the ex- 
ception of the model with — 2.35, which experiences some 
threshold-induced star-formation gasping at times > 12 Gyr. 
The reason why the models with ft = 1 and j3 = 2 are not 
influenced by the threshold is that in these two cases the 
IGIMF is flatter than the standard IMF (Fig. 1), which im- 
plies larger mass ejection rates and consequently larger gas 
masses available for star formation. The type la and type II 
SN rates computed with f3 = 1 and /3 = 2 are higher than 
those computed with the standard IMF because the SFR 
values are slightly larger at any time. On the other hand, in 
the case with /3 = 2.35 the SN rates are lower than in the 
standard case because the IGIMF is slightly steeper than 
the standard IMF in the mass range of the type la and type 
II SN progenitors. 

The evolution of the gas surface mass density is weakly 
sensitive to the adopted value of /3. As explained above, a 
lower j3 implies a higher mass ejection rate from dying stars, 
consequently a larger gas mass at any time. The fact that 
the stellar mass surface density computed with the stan- 
dard IMF flattens at times > 10 Gyr is basically due to the 
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effect of the star formation threshold, which inhibits star 
formation at late evolutionary times. This flattening is not 
predicted for the three cases computed adopting the IGIMF 
since, as we have already seen, the star formation histories 
are not strongly influenced by the SF threshold. 

In Fig. 3, we show the predicted time evolution of the 
[Fe/H] (lower panel) and the [O/Fe] - [Fe/H] plot. In general, 
lower values of P, corresponding to flatter IGIMFs, produce 
at any time higher [Fe/H] values and higher [a/Fe] values 
at a given metallicity. 

4.3 Local observables and their dependence on 
model parameters 

The observables considered in this paper, summarized in 
Tab. 2, are: 

• the abundance ratios, observed in local stars of various 
metallicities; The abundance ratios between two elements 
formed on different timescales are very useful diagnostics, 
since they allow us to constrain the star formation history 
of the studied system. They provide us with information on 
the relative roles of various stellar sources in the chemical 
enrichment of the interstellar medium (Matteucci 2001). In 
particular, the study of [a/Fe]^ is of major importance, ow- 
ing to the difference in the timescales for a-elements and Fe 
production. The abundance ratios are sensitive to the as- 
sumption of the IMF and to the star formation history. 

• the solar abundances, which are useful to test whether 
our models correctly reproduce the metallicity of the Sun at 
the epoch when the Solar System formed, i.e. ~ 4.5 Gyr ago; 
The solar abundances provide information on the integrated 
star formation history. 

• the observed stellar metallicity distribution, which tells 
us the differential distribution of the living stars as a func- 
tion of [Fe/H]; 

• the present-day type la and type II SN rates, which are 
sensitive to the star formation history, to the stellar IMF 
and to the fraction of binary systems able to produce type 
la SNe; 

• the present-day gas surface mass density, sensitive 
mostly to the star formation history 

• the present-day stellar surface density, depending 
mainly on the star formation history 

• the present-day mass function, depending on the com- 
bination of the initial mass function and the star formation 
history. 

In the following, we will present all our results obtained with 
the IGIMF assuming three different values for the index P 
of the star cluster mass function. The results obtained with 
the IGIMF are compared to those obtained by means of the 
standard IMF. In each case, we aim at obtaining the best 
match between our models and the set of observational con- 
straints by varying the star formation efficiency v . In every 

t All the abundances between two different elements X and Y 
are expressed as [X/Y] = log(X/Y) - log(X/Y) Q , where (X/Y) 
and (X/Y)q are the ratios between the mass fractions of X and 
Y in the ISM and in the Sun, respectively. We use the set of solar 
abundances as determined by Grevesse at al. (2007). 



single case for P, the best model is the one providing the 
best match simultaneously to the calculated SN rates, gas 
and stars mass surface densities, solar abundances, abun- 
dance ratios and present-day mass function. 

We do not consider the infall timescale as a free pa- 
rameter. Our assumption is based on the fact that numeri- 
cal dynamical models for disc galaxy formation indicate in- 
fall time-scales of several Gyr (Larson 1976; Samland et al. 
1997). An infall timescale of 5-6 Gyr is indicated also by 
cosmological SPH simulations of disc galaxies in a standard 
A-cold dark matter cosmology (Sommer-Larsen et al. 2004). 
However, we have tested the effects of varying the infall 
timescale, verifying that variations of the e-folding time r 
of 1-2 Gyr have a negligible impact on our results. 

4.4 P = l 

In Fig. 5, we show the star formation history, the time evo- 
lution of the type la and type II SN rates, of the gas surface 
density and of the stellar mass density for three models with 
P — 1 and assuming various SF efficiencies v, compared to 
the results obtained with the standard IMF. By assuming 
P = 1 for the embedded cluster mass function, the model 
which best reproduces our set of observational constraints is 
characterized by a SF efficiency v = 0.1 Gyr -1 . This can be 
seen from Fig. 4, where we show how the "fitness" quantity 
defined in section 4 behaves as a function of p. Among the 
models with 0=1, the one with the highest fitness value is 
that with v = 0.1 Gyr -1 . 

From Fig. 5 it is clear that there are some observables 
which cannot be reproduced accurately, such as the present 
gas and stellar mass densities. Furthermore, All the models 
with P = 1 tend to overestimate the type II SN rate. This 
result should be expected, since the assumption p = 1 pro- 
duces an IMF flatter than the standard one, hence richer in 
massive stars. 

In Table 3, we present the predicted solar abundances 
for the most important heavy elements, obtained with the 
standard IMF and with the best models using the IGIMFs 
for various p. We limit our calculations to the cases of the 
most important chemical elements, for which the standard 
model reproduces the observed abundance pattern with the 
most accurate precision. 

In the case P = 1, the IGIMF being flatter than the 
standard IMF implies heavy element abundances higher 
than the solar ones. By decreasing further the SF efficiency, 
it is possible to improve the match between the observed so- 
lar abundances and the predicted ones, but at the expense 
of the match of the local gas and star surface mass densities. 

In Fig. 6, we present the abundance ratios as a function 
of metallicity, traced by [Fe/H], computed by means of our 
standard chemical evolution model and compared to three 
models with P = 1 and different SF efficiencies. As expected, 
owing to the excess of massive stars and to the very high type 
II SN rates, the predictions obtained with the IGIMF and 
P — 1 overestimate all the [a/Fe] ratios. It is also interest- 
ing to note that the results are weakly dependent on the SF 
efficiencies for values in the range 0.1 < ^/Gyr -1 < 0.5. 
In Fig. 7, we show instead the stellar metallicity distribu- 
tion computed by means of our standard IMF and with the 
IGIMF, showing results for three different SF efficiencies. 
The standard IMF allows us to reproduce the observed stel- 
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Figure 2. From top-left, anti-clockwise: time evolution of the type la SNR, type II SN rate, SFR, surface density in stars and in gas 
computed with the standard IMF (solid lines) and with the IGIMF assuming j3 = 1 (dotted lines), /3 = 2 (dashed lines) and = 2.35 
(dot-dashed lines). For all the models we have assumed the same value for the star formation efficiency, v = 1 Gyr -1 . 




Figure 3. Time evolution of the interstellar [Fe/H] (lower panel) and [0/Fc]-[Fe/H] relation (upper panel) computed with the standard 
IMF (solid lines) and with the IGIMF assuming /3 = 1 (dotted lines), /3 = 2 (dashed lines) and /3 = 2.35 (dot-dashed lines). For all the 
models we have assumed the same value for the star formation efficiency, v = 1 Gyr -1 . 



© — RAS, MNRAS 000, 1-22 



Effects of the integrated galactic IMF on the chemical evolution of the solar neighbourhood 9 



Table 3. Fractional mass in the Sun for various elements as observed by Grevesse et al. (2007, second column) and as predicted by 
means of our model by assuming the standard IMF (third column) and for our best models adopting the IGIMF computed with = 1 
(fourth column), = 2 (fifth column), and = 2.35 (sixth column). 



Element 


Obs. 


standard 1Mb 


1G1M1 


IGIMI 


1G1M1 




(G07) 




= 1 


= 2 


= 2.35 


H 


0.7395 


0.735 


0.732 


0.721 


0.730 


He 


0.2485 


0.253 


0.251 


0.260 


0.259 


C 


2.18(-3) 


1.83(-3) 


1.72(-3) 


2.57(-3) 


2.47(-3) 


O 


5.41(-3) 


5.5(-3) 


8.819(-3) 


8.20(-3) 


3.35(-3) 


Mg 


6.01(-4) 


7.0(-4) 


8.89(-4) 


1.10(-3) 


4.94(-4) 


Si 


6. 70 (-4) 


7.6(-4) 


9.19(-4) 


1.20(-3) 


6.35(-4) 


Fe 


1.17(-3) 


1.21(-3) 


1.19(-3) 


1.96(-3) 


1.32(-3) 


Z 


0.012 


0.0126 


0.016 


0.019 


0.010 



0.65 



0.6 

ro 
m 

C 

0.55 



0.5 



1 1.5 2 2.5 

Figure 4. Open circles: fitness (defined in Section 4) as a function of for various models, each one characterized by a particular star 
formation efficiency v, indicated by the number beside each open circle. The horizontal line indicates the fitness value computed for 
the standard model. From this figure it is clear that the model with = 2 and v = 0.5 is the one providing the best fit to the set of 
observables studied in this paper. 

lar metallicity distribution (SMD) with good accuracy, con- 
cernig either the low-metallicity tail and the peak metal- 
licity value. In Fig. 7, we show both the weighted and re- 
constructed SMDs as observed by Jorgensen (2000). The 
error bars are provided only for the reconstructed SMD and 
are also plotted in Fig. 7. By means of our best model 
computed with the IGIMF and /3 = 1, the predicted SMD 
is not well reproduced. The peak metallicity is underesti- 
mated, and the number of stars with metallicity in the range 
-0.5<[Fe/H]< —0.25 is overestimated. On the other hand, 
the high-metallicty tail is remarkably in disagreement with 



o o.i 

8 8:1s 



O 0.5 
O 0.3 



O 1 - 
O 1.5 



the observations. Models characterised by higher star for- 
mation efficiencies lead to an overabundance of stars with 
metallicity [Fe/H]>0 with respect to the observations. 
We can conclude that, by assuming (3 = 1, no model can 
satisfactorily reproduce at the same time all the observables 
considered in this work. 

4.5 $ = 2 

In Fig. 8, we present the star formation history, the time 
evolution of the type la and type II SN rates, of the gas sur- 
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Figure 5. From top-left, clockwise: calculated time evolution of the star formation history, of the gas surface density, stellar surface 
density and SN rates computed by means of the solar neighbourhood model with the standard IMF (solid lines) and by means of three 
models with the IGIMF, in the case (3 = 1 and assuming three different SF efficiencies: v = 0.1 Gyr -1 (dotted lines), v = 0.25 Gyr -1 
(dashed lines), and v = 0.5 Gyr -1 (dash-dotted lines). In the panel showing the SN rate evolution, the lower and upper curves represent 
the calculated type la and type II SN rates, respectively. Observational data are reported in Tab. 2. 



face density and of the stellar mass density for three models 
with /3 = 2 and different SF efficiencies u, compared to the 
results with the standard IMF and to the available local ob- 
servations. The case with /3 = 2 is more promising than the 
previous one for reproducing the constraints considered in 
this work. In this case, the model which best reproduces our 
set of observational constraints is characterized by a SF ef- 
ficiency v = 0.5 Gyr -1 (see Fig. 4), althogh also the model 
with SF efficiency v — 0.3 Gyr -1 provides satisfactory re- 
sults, as shown by the fitness test. We have also tested mod- 
els with SF efficiency v > 0.5 Gyr -1 , finding fitness values 
lower than the one found with v = 0.3 — 0.5 Gyr - . 

Both type la and II SNRs are larger than those pre- 
dicted with the standard IMF. However, both the present- 
day SN rates predicted with the IGIMF are consistent with 
the local observations by Cappellaro (1996). The gas and 
stellar surface densities predicted with the best model with 
the IGIMF are in agreement with the observed values. 

In the fourth column of Table 3, we present the solar 
abundances predicted with the best model with the IGIMF 
and P = 2. As can be seen from Fig. 8 and Table 3, a satis- 
factory match between the predictions and the observations 
is achievable in this case, although the heavy element abun- 
dances are slightly overestimated. 

In Fig. 9 we show the abundance ratios as a function of 
[Fe/H], compared to available observations in local stars. In 
this case, the main trends of the observed abundance ratios 
as a function of [Fe/H] are well accounted for, as well as the 
abundance ratios at the solar [Fe/H]. 



The results obtained with this choice of /3 are quite sim- 
ilar to the ones achieved with the standard IMF. This result 
could be expected since, as shown in Fig. 1, in the interme- 
diate case with /3 — 2 the IGIMF is very similar to the stan- 
dard IMF. However, the IGIMF is slightly more top-heavy 
than the standard IMF. For the best model, this translates 
into [a/Fe] ratios sligthly larger than those computed with 
the standad IMF. 

As shown in Fig. 10, it is possible to obtain good re- 
sults also for the SMD, with very little difference between 
the model characterized by v = 0.3 Gyr -1 and /3 = 2 and 
the standard model. On the other hand, the model with the 
IGIMF v — 0.5 Gyr -1 is slightly overestimating the number 
of high metalicity stars. However, globally, with ft = 2 sat- 
isfactory results can be achieved also concerning the study 
of the stellar metallicity distribution. 

We can conclude that the assumption of ft = 2 allows 
us to satisfactorily reproduce the set of observational con- 
straints considered in this work. This is an important result, 
given the fact that the IGIMF is computed from first prin- 
ciples. Furthermore, this result allows us to constrain the 
embedded cluster mass function. In the case this function 
may be approximated by a power law, here we have shown 
that the index /3 = 2 is to be favoured with respect to the 
case ;8 = 1. 

The results with /3 = 2 are very similar to those ob- 
tained with the standard IMF, and in principle, on the basis 
of the results described in this section, it may be difficut to 



© — RAS, MNRAS 000, 1-22 



Effects of the integrated galactic IMF on the chemical evolution of the solar neighbourhood 11 




-3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 

[Fe/H] 




Figure 6. Predicted [0/Fc]-Fe/H] (upper panel) and [Si/Fe]-Fe/H] (lower panel) computed by means of the solar neighbourhood model 
with the standard IMF (solid lines) and by means of three models with the IGIMF, in the case /3 = 1 and assuming three different SF 
efficiencies: v = 0.1 Gyr - 1 (dotted lines), v = 0.25 Gyr -1 (dashed lines), and v = 0.5 Gyr -1 (dash-dotted lines). The predictions are 
compared to observational data from various authors. Asterisks: Cayrel et al. (2004); open circles: Bensby et al. (2003); solid diamonds: 
Francois et al. (2004). The solid squares are from a compilation of data by Cescutti (2008). 



discriminate between the two scenarios. In Sect. 4.8, we will 
suggest the use of a diagnostic allowing us to put further 
constraints on the IMF in the solar neighbourhood and to 
disentangle between the standard IMF and the IGIMF. 

4.6 ft = 2.35 

The assumption of ft = 2.35 produces fewer massive stars 
than in the standard case, and this has some impact on the 
predicted SN rates and gas and stellar surface mass densi- 
ties. In Fig.s 11, 12, and 13 we show the predicted time 
evolution of the SFR, SN rates, gas and stellar mass densi- 
ties computed by assuming the IGIMF with ft = 2.35 and 
SF efficiencies v — 1 Gyr" 1 , v — 1.5 Gyr -1 , and v = 2 
Gyr -1 , respectively, compared to the standard model and 
to the local observational values. In the sixth column of Ta- 
ble 3 we show the solar abundances predicted with the best 
model with ft = 2.35, characterized by a SF efficiency v = 2 
Gyr -1 (see Fig. 4). This implies a quicker gas consumption 
timescale and stronger effects of the SF threshold than in the 
standard case, which has a SF efficiency v — 1 Gyr -1 . While 
the model with the standard IMF is influenced by the SF 
threshold only at times > 11 Gyr, the adoption of a higher 
SF efficiency causes strong threshold effects already after 1 
Gyr of evolution. This is visible mainly in the predicted SF 
history, type II SNR and gas surface density evolution. 

For some elemnents (0,Mg), the model reproducing 
best the data of Fig. 11 provides solar abundances lower 
than the observed ones (see Tab. 3). On the other hand, as 



shown by Fig. 14, the analysis of the abundance ratios as 
a function of metallicity indicates [a/Fe] values lower than 
the ones computed with the standard IMF. Finally, in Fig. 
15, we show the SMD computed with the IGIMF assuming 
ft = 2.35 compared to the observations and the results of the 
standard IMF. In this case, the peak metallicity computed 
with the IGIMF is shifted leftwards by 0.2 dex with respect 
to the observations and to the standard model. From the 
results discussed in this section we can conclude that by as- 
suming ft = 2.35 it is not possible to reproduce at the same 
time all the observational constraints considered in this pa- 
per. 

4.7 Effects of other parameters: the star 
formation threshold and Ai a 

At this stage, it may be interesting to study the effects which 
other important parameters have on our results, such as the 
SF threshold and the type la SN realization probability Ai a . 

As shown in Sect. 4.5, the models which best repro- 
duce the constraints considered in this paper are not sensi- 
tive to the adopted SF threshold. This is visible from the 
smooth behaviour of the calculated SF history, whereas we 
have seen that when the effect of the threshold is important, 
it presents an oscillatory behaviour, as happens for the SFH 
of the model with the standard IMF. The case with ft = 1 is 
characterized by a high amount of mass restored by massive 
stars into the ISM, keeping its density always significantly 
above the threshold. Furthermore, in this case the SF effi- 
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Figure 7. Present-day stellar metallicity distribution in the Solar Neighbourhood. The solid line, the dotted line, the dashed line and 
the dash-dotted line represent the predicted stellar metallicity distribution computed with the standard IMF and with the IGIMF in 
the case /? = 1 and SF efficiency v = 0.1 Gyr - 1 , v = 0.25 Gyr -1 , and v = 0.5 Gyr — , respectively. The thin and thick solid histograms 
represent the observational weighted and reconstructed SMD from Jorgensen (2000), respectively. For the reconstructed distribution, 
error bars are plotted. 



ciency which allows us to best reproduce the local features 
is lower than in the standard case. This causes the gas mass 
to be less efficiently consumed, so that the gas density level 
is always considerably above the threshold. The same is true 
for the case with f3 = 2. In these two cases, halving the SF 
threshold to 4Mq/pc 2 does not produce any effect on our 
results. On the other hand, the adoption of a significantly 
higher threshold, say > 10M© /pc 2 would be at variance with 
the observations in local spiral discs (Martin & Kennicutt 
2001). The case with f} — 2.35 is the most sensitive to the as- 
sumption of the threshold. We show the effects of a threshold 
4Mq/pc 2 on the star formation history in Fig. 16. The de- 
crease of the threshold causes the SF history to be smoother 
and to enhance gas consumption, therefore at the present 
time the galaxy model with EtAr = 4M© /pc 2 presents a gas 
density and a stellar density slightly lower and higher than 
the model with E t h r = 7Mq/pc 2 , respectively. On the other 
hand, the assumption of a lower threshold has no effect on 
the abundance ratios, and this can be seen in Fig. 17. 

It is worth noting that recent GALEX results indicate 
that the SF cutoff visible from the Ha profiles of local star 
forming galaxies is instead absent in UV profiles (Boissier 
et al. 2007), suggesting that the existence of a SF threshold 
may be an observational selection effect, naturally explained 
by the IGIMF theory as shown by Pflamm-Altenburg & 
Kroupa (2008). However, as shown by chemical evolution 
results, the SF threshold is fundamental in reproducing the 
metallicity gradients observed in the MW and in local galax- 
ies (Chiappini et al. 2001; 2003), unless a variable star for- 



mation efficiency through the disc is assumed (Colavitti et 
al. 2009). 



It may be also interesting to discuss the possible effects 
of the type la SN realization probability Ai a on our results. 
A change in Ai a has effects on the predicted type la SN rate, 
on the [X/Fe]-[Fe/H] plots and on the SMD plot, whereas 
it leaves unchanged any other quantity studied here. In the 
case j3 = 1, the model which best reproduces our set of ob- 
servational constraints is characterized by a SF efficiency 
v = 0.1 Gyr - . To reproduce the abundance ratios as a 
function of [Fe/H], a higher value of Ai a would be required, 
so that all the curves would move downwards. This would 
improve also the fit to the observed stellar metallicity distri- 
bution, since an increase of the type la SN rate efficiciency 
would shift the computed SMD rightwards, probably in bet- 
ter agreement with the observations. However, as shown in 
Fig. 5, even with a higher value for Ai a , the observed local 
type II SN rate, the stellar mass density and the gas density 
(which do not depend on Ai a ) would still not be reproduced. 
We have seen that in the case f} — 2, no substantial modifi- 
cation of Aia is required. 

Finally, in the case /3 = 2.35, a lower value of Ai a would 
be required to improve our match to the observations in the 
[X/Fe]-[Fe/H] plot, however this would shift the predicted 
SMD leftwards, exacerbating the discrepancy between the 
computed SMDs and the the observed distribution. 
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Figure 8. From top-left, clockwise: predicted time evolution of the star formation history, of the gas surface density, stellar surface 
density and SN rates computed by means of the solar neighbourhhod model with the standard IMF (solid lines) and by means of three 
models with the IGIMF, in the case fj = 2 and assuming three different SF efficiencies: v = 0.1 Gyr -1 (dotted lines), v = 0.3 Gyr -1 
(dashed lines), and v = 0.5 Gyr -1 (dash-dotted lines). Symbols as in Fig. 5. 



4.8 The Present-day Mass Function 

The present-day mass function (hereinafter PDMF) repre- 
sents the mass function of living stars as observed in the so- 
lar neighbourhood. This quantity is an important diagnostic 
since it allows us to test, beside the IMF, the star formation 
history of our best model, providing further information on 
our parameters complementary to the ones previously dis- 
cussed. 

The number of stars formed in the time interval (t, t+dt) 
and in the mass interval (m, m + dm) is 

£(m,t)i>(t)dmdt, (13) 

(Tinsley 1980). The present number per unit mass of stars 
with lifetime r m < To is 

/•To 

n(m)= / i>(t)£(m,t)dt. (14) 

JT -r{ m ) 

The total number of stars with mass between M m i n = 
0.1 Mq and m is N(m), calculated as 

N(m) = I I ip(t)^(m',t)dtdm, (15) 

''Mmin Jt — r(m') 

where To is the present time, equal to 14 Gyr. The present- 
day mass function can be calculated as 

PDMF — dN(m)/d(logm). (16) 

In Fig. 18, we show the local PDMF as estimated observa- 
tionally by various authors, and as predicted by means of our 
models assuming the standard IMF and the IGIMF for all 



the three different cases for /3, assuming the SF efficiencies 
of the best models as indicated by our fitness test described 
in Sect. 4.1. The analytical PDMF of Chabrier (2003) has a 
log-normal profile for masses m < IMq: 



dN 



zxp{- 



(log 



log m c y 



d (log m) 

with m c = O.O79M0 and a 
it is given by a power-law: 
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0.69, and for masses m > 1M ( 

if < log(m/M e ) < 0.54, 
if 0.54 < log(m/M e ) < 1.26, 
if 1.26 < log(m/M Q ) < 1.80. 



o 




The analytical PDMF of Kroupa et al. (1993) is a power-law: 

if log(m/M Q ) < -0.3, 
if - 0.3 < log(m/M Q ) < 0, 
if < log(m/M Q ), 

The Kroupa et al. (1993) and Chabrier (2003) PDMF are 
very similar at masses M < IMq, and in the same mass 
range they are in very good agreement with the Miller & 
Scalo (1979) PDMF. For masses M > IMq, no recent up- 
date exists and the reference measures are those of Miller & 
Scalo (1979) and Scalo (1986), which are in very good agree- 
ment in this mass range. In the light of this, for purposes 
of clarity and to avoid confusion, as observational data in 
Fig. 18 we plot only the PDMF as determined by Miller & 
Scalo (1979). This is also the one used for the fitness test 
described in Sect. 4.1. 

From the analysis of Fig. 18, interesting information 
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Figure 9. Predicted [0/Fc]-[Fe/H] (upper panel) and [Si/Fe]-[Fe/H] (lower panel) computed by means of the solar neighbourhood model 
with the standard IMF (solid lines) and by means of three models with the IGIMF, in the case j3 = 2 and assuming three different 
SF efficiencies: v = 0.1 Gyr -1 (dotted lines), v = 0.3 Gyr -1 (dashed lines), and v = 0.5 Gyr -1 (dash-dotted lines), compared to 
observational data from various authors (see caption of Fig. 6). 
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Figure 10. Present-day stellar mctallicity distribution in the Solar Neighbourhood. The solid line, the dotted line, the dashed line and 
the dash-dotted line represent the predicted SMDs computed with the standard IMF, with the IGIMF in the case j3 = 2 and SF efficiency 



v = 0.1 Gyr~ 



0.3 Gyr 1 , and u = 0.5 Gyr 1 , respectively. The histograms are the observational data described in Fig. 7. 
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Figure 11. From top-left, clockwise: predicted time evolution of the star formation history, of the gas surface density, stellar surface 
density and SN rates computed by means of the solar neighbourhhod model with the standard IMF (solid lines) and by means of the 
best model with the IGIMF, in the case /3 = 2.35 and a SF efficiency v = 1 Gyr -1 (dotted lines). In the panel showing the SN rate 
evolution, the lower and upper curves represent the predicted type la and the predicted type II SN rates, respectively. Symbols as in 
Fig. 5. 



can be drawn on the shape of the IMF and on the effects 
of other parameters such as the SF threshold. The PDMF 
computed with the standard IMF is in excellent agreement 
with the observational data in the range 0.4 Mq - 2 Mq. At 
stellar masses lower than 0.4 Mq, the standard IMF is too 
steep and overestimates the observed number of low mass 
stars. On the other hand, the model computed with the stan- 
dard IMF underestimates the observed distribution of mas- 
sive stars with masses > 30Mq . This is due to the adoption 
of the SF threshold and its strong effects on the SF history of 
the solar neighbourhood at late times, inhibiting recent SF 
and hence causing the underabundance or absence of very 
massive stars. 

The best models calculated with the IGIMF for /3 = 2 
and P = 2.35 have different values for the SF efficiencies, but 
provide all similarly a very good fit to the observed PDMF. 
The model with /3 = 1 slightly underestimates the low-mass 
end of the PDMF. 

The small offsets among the computed distributions are due 
to different values adopted for the star formation efficiencies. 
These assumptions produce different values for the present- 
time stellar mass densities, which reflect into the small off- 
sets visible in Fig. 18. 

Our results indicate that the values /3 = 2 and /3 = 2.35 
should be preferred. In the case P = 2.35, apparently the SF 
threshold does not affect the computed PDMF. The best 
model with the IGIMF and /3 = 2.35 presents no truncation 
in the PDMF since, even if the SFR has a gasping behaviour, 
the frequency of recent star formation events is sufficient to 



guarantee a significant presence of high-mass stars still liv- 
ing at the present time. 

Finally, it is important to stress the importance of the fitness 
test used to determine which is the best model to reproduce 
all the obsevables, and how its results depend on the list 
of observables taken into account. In Table 4, we present 
our results for the fitness test performed considering (i) the 
whole set of observables and, for comparison, (ii) computed 
by considering all observables but the PDMF. In column 1, 
for each model we report the j3 values. In column 2, we show 
the SF efficiencies. In column 3, we show the fitness obtained 
by using all the obsevables, and in column 4 we present the 
fitness computed by excluding the PDMF. It is clear that, 
once the PDMF is excluded and, if one considers the abil- 
ity to reproduce all the chemical evolution constraints and 
physical properties such as the SN rates or the gas mass den- 
sity, the two models with ft = 2, v = 0.3 Gyr -1 and v = 0.5 
Gyr" 1 are the best, but the standard model provides com- 
parable results. However, the uncapability of the standard 
model of accurately reproducing the observed PDMF has a 
strong impact on the fitness computed considering the whole 
set of observables. In this case, the two models with /3 = 2 
yield better results than the standard model. 

4.9 The time variation of the IGIMF 

In Fig. 19, we show how the IGIMF varies as a function of 
time in the case of the three best models, with /3 = 2.35 
and v = 2 Gyr -1 (upper panel), with /3 = 2 and v = 0.5 
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Figure 12. From top-left, clockwise: predicted time evolution of the star formation history, of the gas surface density, stellar surface 
density and SN rates computed by means of the solar neighbourhhod model with the standard IMF (solid lines) and by means of the 
best model with the IGIMF, in the case /3 = 2.35 and a SF efficiency v = 1.5 Gyr -1 (dotted lines) In the panel showing the SN rate 
evolution, the lower and upper curves represent the predicted type la and the predicted type II SN rates, respectively. Symbols as in 
Fig. 5. 



Gyr -1 (central panel) and with /3 = 1 and v — 0.1 Gyr -1 
(lower panel). In each panel, we show the IGIMF at various 
epochs, each one characterised by different SFR values. 

The best model with /3 — 1 is the one characterized by 
the lowest SF efficiency and presents a rather smooth SFH, 
although it has stronger variations than the other cases. In 
fact, as can be seen from Fig. 5, at time 1 Gyr the model 
with v = 0.1 Gyr -1 has a SFR < lM@/yr, whereas at the 
present time it has a SFR ~ 2Ms,/yr. On the other hand, in 
the period going from 1 Gyr to the present time, the other 
two models have SFR variations within a factor of 2, smaller 
than in the case with u — 0.1 Gyr -1 . The variations in the 
SFH of the best model are the cause of the time variations 
in the IGIMF visible in Fig. 19 at the highest stellar masses 
(log(M/M ~ 1.4)). 

The best model with /3 = 2 presents a rather smooth 
SFR characterised by no strong variation in the period 1 
Gyr-14 Gyr (see Fig. 8), not influenced by the effects of the 
star formation threshold. This reflects the very small varia- 
tion of the IGIMF with cosmic time. 

The threshold plays an important role in the evolution of 
the model with f} — 2.35, whose SFH presents a gasping 
behaviour throughout the whole cosmic time. As a conse- 
quence, the IGIMF presents strong variations. 

Our results show that within the IGIMF theory, if the 
absolute best model reproducing the properties of the solar 
neighbourhood is the one with /3 = 2 and v = 0.5 Gyr -1 , 
the variation of the IMF with time must be small. 

In principle, other objects with extreme variations in 



the SFH such as starburst galaxies should exhibit a strongly 
time-varying IMF. This will be an interesting subject for 
future work. 



5 CONCLUSIONS 

The initial mass function regulates the number of stars of 
different intial mass born per stellar generation, hence it 
plays a fundamental role in galactic chemical evolution stud- 
ies. The aim of this paper was to test the effects of adopting 
the integrated galactic initial mass function (IGIMF) on the 
chemical evolution of the solar neighbourhood. The IGIMF 
is computed from the combination of the stellar intial mass 
function, i.e. the mass function in single star clusters, and 
the embedded cluster mass function, i.e. a power law with 
index /3, and taking into account that within each single 
cluster, the maximum stellar mass is a function of the total 
mass of the cluster. The result is a time- varying IMF which 
is a function of the galactic star formation rate. We applied 
the formalism developed by Weidner & Kroupa (2005) to a 
chemical evolution model for the solar neighbourhood, based 
on the two-infall model by Chiappini et al. (1997). For the 
embedded cluster mass function, we tested three different 
values of j3 and studied various physical quantities and abun- 
dances for various chemical elements, comparing our results 
to the ones obtained by means of a standard IMF, constant 
in time and similar to the Scalo (1986) IMF. A statistical 
test to determine which is the best model in reproducing the 
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Figure 13. From top-left, clockwise: predicted time evolution of the star formation history, of the gas surface density, stellar surface 
density and SN rates computed by means of the solar neighbourhhod model with the standard IMF (solid lines) and by means of the 
best model with the IGIMF, in the case = 2.35 and a SF efficiency v = 2 Gyr -1 (dotted lines). In the panel showing the SN rate 
evolution, the lower and upper curves represent the calculated type la and type II SN rates, respectively. Symbols as in Fig. 5. 



Table 4. Fitness test for various models in two different cases. In columns 1 and 2, for each model the value of and SF efficiency values 
are reported, respectively. In columns 3 and 4, the results of the fitness test computed by considering the whole set of observables and 
all observables but the PDMF are shown, respectively. 



v Fitness, Fitness 

all observables all obs. except PDMF 



1.00 


0.10 


0.54 


0.56 


1.00 


0.25 


0.53 


0.57 


1.00 


0.50 


0.53 


0.55 


2.00 


0.10 


0.47 


0.51 


2.00 


0.30 


0.66 


0.69 


2.00 


0.50 


0.68 


0.72 


2.35 


1.00 


0.48 


0.53 


2.35 


1.50 


0.46 


0.55 


2.35 


2.00 


0.52 


0.55 


Standard 


1.00 


0.56 


0.68 



set of observational data considered in this work has been 
developed. Our results can be summarized as follows: 
1) The value of f3 has important effects on the predicted 
star formation history. Also the effects of the SF threshold 
may be weaker or stronger, depending on the value chosen 
for /?. In general, lower absolute values of /3 imply a flatter 
IGIMF, hence a larger number of massive stars and larger 
mass ejection rates. This translates into a larger amount of 
gas available for star formation and gas density values never 



below the threshold and a smooth star fomation history. 

2) The value of P has an obvious deep impact on the pre- 
dicted SN rates. Beside this, also other quantities can be 
influenced by this parameter. In general, a lower f3 implies 
higher mass ejection rates from massive stars, hence at any 
time a larger gas mass density. 

3) The interstellar abundances are strongly influenced by 
the parameter /3. In general, lower absolute values of ft im- 
ply higher metallicities and higher [O/Fe] values at a given 
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Figure 14. Predicted [0/Fc]-[Fc/H] (upper panel) and [Si/Fc]-[Fe/H] (lower panel) computed by means of the solar neighbourhood 
model with the standard IMF (solid lines) and by means of three models with the IGIMF, in the case /3 = 2.35 and assuming three 
different SF efficiencies: v = 1 Gyr — 1 (dotted lines), v = 1.5 Gyr — 1 (dashed lines), and v = 2 Gyr -1 (dash-dotted lines), compared to 
observational data from various authors (see caption of Fig. 6). 



metallicity. 

4) We have considered several chemical evolution con- 
straints, including the observed local SN rates, local gas and 
stellar surface densities, the abundance ratios in local stars 
and the stellar metallicity distribution and, by varying the 
SF efficiency, we tested which assumption for the embedded 
cluster mass function exponent j3 provides the best simul- 
taneous match to all of these observables. Our fitness test 
indicates that the best result has been obtained by assum- 
ing j3 = 2, which produces an IMF resembling that of our 
standard model and that allows us to account for most of 
the observables. 

5) Useful hints on the initial mass function come from the 
study of the present day mass function, i.e. the mass function 
of living stars observed today in the solar neighbourhood. 
The PDMF is a very good test of the chosen star formation 
history (Elmegreen & Scalo 2006) once an IMF has been 
assumed and a valuable test for the IMF once the SFR is 
assumed. In the case of the standard IMF, the star fomation 
history assumed here predicts a present-time SFR in good 
agreement with the observations, and is based on the results 
by Kennicutt (1989). Models with the IGIMF and different 
assumptions for /3 provide different results. Lower values for 
/3 produce a PDMF richer in massive stars, whereas higher 
values of /? imply steeper IGIMFs and a lower relative num- 
ber of massive stars. The model with /3 = 2 should be con- 
sidered the best since it allows to reproduce at the same 
time the observed PDMF and to account for most of the 
chemical evolution constraints considered in this work. 

In the future, we plan to extend our study of the effects 



of the IGIMF in galaxies of different morphological types. As 
forthcoming step, we will investigate how the IGIMF affects 
the chemical evolution of dwarf galaxies. 
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Figure 15. Present-day stellar metallicity distribution in the Solar Neighbourhood. The solid line, the dotted line, the dashed line 
and the dash-dotted line represent the predicted SMDs computed with the standard IMF, with the IGIMF (in the case = 2.35) and 
SF efficiency u = 1 Gyr" 1 , with the IGIMF and SF efficiency v = 1.5 Gyr" 1 , and with the IGIMF and SF efficiency v = 2 Gyr" 1 , 
respectively. The histograms are the observational data described in Fig. 7. 
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Figure 16. From top-left, clockwise: predicted time evolution of the star formation history, of the gas surface density, stellar surface 
density and SN rates computed by means of the model with j3 = 2.35 and a SF efficiency v = 1.5 and adopting a SF threshold 



^thr = 7M@pc (thin dashed lines) and S t 



AMqpc (thick dot-dashed lines). 
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Figure 17. Predicted [0/Fe]-[Fc/H] (upper panel) and [Si/Fc]-[Fc/H] (lower panel) computed by means of the model with /3 = 2.35 and 
a SF efficiency v = 1.5 and adopting a SF threshold £th r = 7Mqpc~ 2 (dashed lines) and £thr = 4M0pc -2 (thick dot-dashed lines), 
compared to observational data from various authors. 
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Figure 18. Predicted present-day mass function computed by assuming the standard IMF (thick solid line) and the IGIMF with = 1 
(thick dotted line), (3 = 2 (thick dashed line) and j3 = 2.35 (thick dot-dashed line). The observational PDMF is the one of Miller & Scalo 
(1979) (solid triangles with error bars). 
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Figure 19. Evolution of the IGIMF as a function of the cosmic time for the models with /3 = 2.35 and u = 1.5 Gyr -1 (upper panel), 
with /3 = 2 and v = 0.5 Gyr -1 (central panel) and for the model with /3 = 1 and u = 0.25 Gyr -1 (lower panel). In each panel, the curves 
represent IGIMFs at various times, calculated on the basis of the time variation of the star formation history of each model. 
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